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STATISTICAL DISCRETE-SOURCE MODEL OF LOCAL COSMIC RAYS 


R. Ramaty, D.V» Reames, NASA, Goddard Space Flight Cntr., Greenbelt, Md. 

R.E. Lingenfelter*, UniT. of California, Los Angeles 

We wish to investigate the implications of a new upper limit on the sid- 
eral anisotropy of the cosmic rays. After correcting for the influence of the 

interplanetary magnetic field, Elliot et al.^^^ find that cosmic rays of ener- 
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gies ~10 ev have an anisotropy less than 2 x 10 . This value seems to con- 
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flict with estimates of the cosmic-ray lifetime of 1 to 2 x 10 years, 
since for such an anisotropy, the time required for escape from the disk is 
about 3 X 10^ years. 

(3) 

It has been suggested that this discrepancy could be due to fluctuations 
resulting from the statistical nature of the cosmic-ray source distribution 
which at the present just happens to give an extremely low anisotropy at the 
earth. However, from a detailed statistical study, in which we have treated 
the galactic cosmic ray sources as random events in both space and time, we 
find that within reasonable probabilistic limits, the discrepancy between 
isotropy and escape lifetime of 1 to 2 x 10^ years cannot be accounted for 
sclely by fluctuations. 

Within this statistical discrete-source model, we have considered the 
alternate possibility that the streaming velocity is indeed as low as suggest- 
ed by the observed lack of anisotropy. In this case, because of the long 
residence time of the cosmic rays in the galaxy, the lifetimes of the various 
nuclei are determined principally by nuclear interactions and not by escape. 

In particular, the proton lifetime is much longer than the lifetimes of the 
heavier nuclei. The apparent discrepancy between isotropy and lifetime can 
then be readily understood, since the anisotropy measurements refer to the 
total cosmic ray flux which consists mainly of protons, while the lifetime 
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estimates are based principally on the amount of matter traversed in inter- 
stellar space by nuclei in the meditim group. In this respect, the statistical 
model does not greatly differ from the steady state models which were used in 
previous calculations. However, when the combined effects of the statistical 
model and a low streaming velocity are simultaneously considered, several of 
the measurable properties of the cosmic rays, and in particular the abundances 
of nuclei heavier than iron, will exhibit large fluctuation which must be taken 
into account in order to interpret the observational data. 

In this letter, we wish to present the outline and the results of calcu- 
lations of the fluxes, interactions, lifetimes, and anisotropies of the various 
chemical components of the cosmic rays in a statistical discrete-source model. 
We shall present the details of these computations in a separate paper. 

We have assumed that the cosmic ray sources are supernova explosions con- 
fined to a disk of thickness 200 pc with a random spatial and temporal distrib- 
ution normalized to a local rate which corresponds to 1 supernova per 100 years 
in the galaxy. We treat each supernova as a point source in space and time and 
we assume for simplicity that cosmic ray propagation in interstellar space can 
be characterized by 3-dimensional isotropic diffusion with a constant mean free 

path X. We take into account escape from the disk by introducing a character- 
2 

istic time, * 3d /(2\c), where d ^ 500 pc is the distance in the plane 

along which cosmic rays are conveyed to the surface of the disk by the random 

(3) 

walk of the field lines , and escape from the galaxy. We also take into 

account nuclear disintegration by introducing a destruction lifetime, 

where n^ is the hydrogen density in interstellar space and 

-3 

is the destruction cross section. We have used n^ = 1 cm , and = 30 mb 

2/3 

for protons, and ® 50 A mb for the heavier nuclei. We limit our dis- 
cussion to relativistic nuclei, so that ionization losses are negligible. 
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The energy density at the earth resulting from a supernova at a distance 
r^^ and of age t^, in the form of a particular nuclear component whose des- 
truction lifetime is is given by 

where is the total energy released by the supernova in particles of the kind 
considered. We assume that the values of areihe same for all supernovae and 
we denote the total cosmic ray energy from one supernova by W , 

The energy density, lifetime and anisotropy of one nuclear component from 
all supernovae (i = l,o.o,N) are given by 

UX zz Z (2) 

t z ^ ^ ^ O') 

J r 3/z iir-' i" (I K-v y./ct. fy ( I 

where x^, y^, are the coordinates (x^ith the earth at the origin) of the i-th 
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supernova. We have numerically evaluated equations (1) through (4) for H’^,0 , 
Fe^ , and U nuclei and for a range of the scattering mean free path Using 

a statistical treatment which we shall describe in (^ tall in a subsequent paper, 
we have evaluated the mediatis, m, and the la and^^^le^l s of all quantities of 
interest. These levels for the proton anisotropy are plotted in Figure 1 as a 

^l) „4 

function of X. As can be seen, the upper limit of 2 x 10 , can be best 

xinderstood in terms of a low scattering mean free path, since, if X = 1 pc 
corresponding to ^ 10 years (for d = 500 pc), the probability is very low 

(<o02) that S< 2 X 10"^. In view of the fact that the anisotropy measurements 
only give an upper limit, it is most probable that X ^ 0.1 pc. This value is 
about an order of magnitude smaller than that deduced from studies of wave- 
particle interactions in interstellar space. According to Wentzel (private 
communication), however, these studies did not include the turbulence produced 
by low-energy cosmic rays, which, together with other uncertainties in the 
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properties of the interstellar medium, could substantially decrease the effec- 
tive mean free path. 

The ratio of the light nuclei (3 ^ Z ^ 5) to the medixim nuclei (6 ^ Z ^ 8) 
at relativistic energies is plotted as a function of X in Figure 2. This ratio 
was computed by assuming that the former are produced by the fragmentation of 
the latter in interstellar space. The fragmentation cross section that we have 
- used was 72 mb^^^and the destruction cross sections for the medium and light 
groups were 185 mb and 130 mb^^^respectively. The measured L/M ratio^^^at 
relativistic energies is about 0.25 + 0.05. From Figure 2, we see that within 
the IcT levels of the present model, this ratio is consistent with 10“^ ^X^10’*^pc. 

Also shown in Figure 2 are the medians and the l<j levels of the lifetimes 

of the 0^^ and nuclei, together with the escape lifetime, At 

large X*s the lifetimes are dominated by escape and are the same for all nuclei. 

At lower values of X 1 pc) nuclear disintegrations become inqportant and for 

differing nuclei yield different lifetimes all of which are less than In 

-4 

the limit of very low X's 10 ), where all cosmic rays come from a single 

nearby source, the lifetimes of the various nuclei again approach a common 
value equal to the age of the source. The L/M ratio measures the medium group 
lifetime which, in the range 10“^pc iS X ^ 10" ^pc is determined principally by 
nuclear destruction rather than escape. Moreover, for these values of X, even 
the proton lifetime is dominated by inelastic processes and is much less than 
the escape lifetime. For the same reason, the lifetimes of the very high 
energy electrons (^20 BeV) would be dominated by radiation losses and not by 
escape. The statiijtical properties of the present model have important effects 
on the electron speccrum at these energies and must be included in any inter- 
pretation of the electron measurements in terms of the spectra of these parti- 
cles at the sources. In principle, the cosmic ray lifetimes could also be 
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directly determined from the measured abundances of certain radioactive nuclei 
such as and • At the present time, however, the measurements 

are not accurate enough for a positive determination. 

The proton energy density at the earth, the composition ratios and the 
number of significant sources are shown in Table 1. Although we have consid- 
ered all supernovae in the galaxy, because of the combined effects of nuclear 
disintegrations, escape, and propagation, only a relatively small number of 
sources contribute a major fraction of the cosmi,; ray ilux at the earth at any 

one time. For protons, the number ef significant sources M is a maximum at 
-2 

X X 10 , and for heavier ntic.lei, M is a maximum at higher values of X* The 

maximum number of sources obvious ly leads to minf.mum fluctuations in the cosmic 
ray flux. For both high and low values of X the ages and the number of signifi- 
cant sources approach a common limit for ail nuclear components. At large X*s 
the cosmic rays are produced by a young and distant source, whereas at small X*s 
by an old and close one. in the intermediate range, the number of significant 
sources for different nuclei, like the lifetimes, are different and decrease 
with increasing nuclear mass. It is this property of the source distribution 
that is required to account not only for the apparent discrepancy between iso- 
tropy and lifetime but also for the interactions and abundances of the heavy 
nuc lei. 

A lower limit can also be placed on X 10“^pc) since for X < 10 
Mq, would approach 1, and the path length distribution for rela- 

tivistic 0, Fe and U nuclei become the same, which is inconsistent with observed 
fragmentation ratios of products of medium and iron group nuclei^^^\ as well 
as the detection of nuclei with Z ~ 

We suggest, therefore, that the diffusion mean free path is probably con- 
fined to the range 10“^pc ^ X ^ 10" ^pc. The upper limit is determined by the 
observed isotropy of the cosmic rays; the lower limit is determined by the re- 



- 6 - 


qulrement that different nuclei have different path length distributions, and 
the range Is consistent with the observed L/M ratio. 

For the range 10" ^pc ^ X ^ 10" ^pc the cosmic ray energy per supernova, 

Wgjj, required to account for the observed energy density at the earth of 10"^^ 
erg cm"^. Is between .3 to lx 10^^ erg. This energy Is well within the range 
of the estimated cosmic-ray output of supemovao. 

For the same range of X, the fluctuations In the 0/H, Fe/H and U/H ratios 
become large and make a comparison between the observed chemical composition of 
the cosmic rays and source compositions, based on models of nucleosynthesis, 
extremely difficult. The reason for this Is our lack of knowledge of the dis- 
tances and ages of the cosmic ray sources, although It has recently been 
(13) 

suggested that the pulsars may provide this Information. 

Nonetheless, the fluctuations shown In Table 1 directly reflect the ex- 
pected spatial and temporal variations. At X * 3 x 10 pc, the +20% fluctua- 
tions In w„ give a measure of the time averaged variations and are consistent 
with cosmic-ray produced activities in meteorites Shorter term variations 
on the other hand would be larger and could be detectable In the activities 

In deep-sea sedlments^^^' Much larger time averaged fluctuations, however, 

are expected for the heavier nuclei. Fe could vi-ry by about 1007. and the varia- 
tion In U may be several times larger. Much more extensive measurements, possi- 
bly based on the analysis of very heavy particle tracks In lunar samples (P.B. 
Price, private communication) are required to further Investigate these poten- 
tial fluctuations. 
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Figure 1. Cosmic ray anisotropy as a function of the diffusion mean free path 
X. The probability that for a given x the anisotropy is less than the value 
corresponding to the 2a+, lo+, m, Icr-, and To- levels is 0.98, 0.83, 0.5, 0.17, 
and 0.02, respectively. 

Figure 2. Cosmic ray lifetimes and the ratio of light nuclei to medium nuclei 
as a function of X. 
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Xable 1. The proton energy density and the oxygen-, iron- and uranium- to-hydrogen ratios. The quantities 
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